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1 INTRODUCTION 

Stellar flares are generally regarded as the stellar analogues 
of solar flares and their X-ray emission often only differs 
from their solar counterparts in magnitude and duration. 
Stellar flare peak temperatures and emission measures can 
be orders of magnitude greater than what is seen on the 
Sun. The durations of the longest lived stellar flares signifi- 
cantly exceed the longest durations seen in solar flares. For 
a detailed compar i son be tween solar and stellar flares, see 
lAschwanden et all ((2008) . They also differ in the fact that 
stellar flares are spatially unresolved, whereas this is only 
the case in the smallest of flare events on the Sun. Thus it 
is not possible to observe directly where on its host star a 
flare is located, or whether parts of the flare are eclipsed by 
the host star. While typical durations of solar flares are very 
much less than the solar rotation period, this is not the case 
with most stellar flares. As a result, the likelihood that a 
stellar flare undergoes a rotational eclipse is much greater 
than that for solar flares. It is therefore natural to expect 
that although solar and stellar flares are probably produced 
by similar processes, their observational signatures may dif- 
fer. 



For a review of the physical mechanisms in stellar flares, 



ABSTRACT 

Most stellar flares' soft X-ray lightcurves possess a 'typical' morphology, which con- 
sists of a rapid rise followed by a slow exponential decay. However, a study of 216 
of the brightest flares on 161 pre-main sequence stars, observed during the Chandra 
Orion-Ultradeep Project (COUP), showed that many flare lightcurves depart from 
this typical morphology. While this can be attributed to the superposition of multiple 
typical flares, we explore the possibility that the time- variable eclipsing of flares by 
their host stars may also be an important factor. We assume each flare is contained 
within a single, uniform plasma density magnetic loop and specify the intrinsic vari- 
ation of the flare's emission measure with time. We consider rotational eclipse by the 
star itself, but also by circumstellar discs and flare-associated prominences. Based on 
this simple model, we generate a set of flares similar to those observed in the COUP 
database. Many eclipses simply reduce the flare's maximum emission measure or de- 
cay time. We conclude therefore that eclipses often pass undetected, but usually have 
only a modest influence on the flare emission measure profile and hence the derived 
loop lengths. We show that eclipsing can easily reproduce the observed atypical flare 
morphologies. The number of atypical modelled flare morphologies is however much 
less than that found in the COUP sample. The large number of observed atypical flare 
morphologies, therefore, must be attributed to other processes such as multiple flaring 
loops. 

we refer the reader to iBenz fc Giideil <|2010h . A significant 
departure from a potential field configuration in a coro- 
nal magnetic field corresponds to a large amount of ex- 
cess energy being held in that field. When reconnection 
events occur, the coronal magnetic field geometry is sim- 
plified (i.e. becomes closer to a potential field configuration) 
and the resulting configuration corresponds to a lower en- 
ergy state. In this process, a large amount of the excess en- 
ergy is converted into the non-thermal motions of electrons 
and ions which spiral down magnetic field lines and impact 
the stellar chromosphere. This can be seen at radio wave- 
lengths as the charged particles emit gyrosychrotron radia- 
tion. As the energetic electrons impact the chromosphere, 
they emit non-thermal Bremsstrahlung radiation at hard 
X-ray wavelengths, as they become thermalised by random 
Coulomb intera ctions. This model is known as the "thick- 
targe t model" (|Neupert|[l96i ; lBrownlll97ll ; iLin fc Hudson! 
1976). This causes chromospheric plasma to be heated and 
evaporated into the corona where it is contained within mag- 
netic loop structures. Through a combination of mostly ra- 
diative losses and heat condu ction back to the photospher e, 
the evaporated plasma cools l|Antiochos fc Sturrocki ri978'). 

In this paper, we consider the morphologies of the soft 
X-ray lightcurves of spatially unresolved stellar flares. In the 
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Figure 1. Example flare showing how eclipsing can cause an atypical flare morphology. The flare is assumed to be contained within a 
single magnetic loop with a uniform plasma density. The example emission measure curve is the flare on COUP source 649 given in Fig. 
[3] The images at the top show the position of the single magnetic loop containing the flaring plasma, the geometry of which has been 
determined from the best fit emission measure curve given in Fig. [3] Flare phase 0.0 and 1.0 represent the beginning and end of the flare 
respectively. 



majority of cases, the morphologies of typical flares can be 
broken down into two distinct phases. The first phase con- 
sists of a rapid increase in luminosity due to the heating and 
evaporation of chromospheric plasma. This is followed by 
a slow exponential decay due to cooling. However, a large 
number o f stellar flares show m ore complex atypical mor- 
phologies ijGetman et al.ll2008ah . Among these atypical mor- 
phologies are flares with longer rise phases and no clear peak 
and flares with multiple peaks or dips in their lightcurves. 
The interpretation of these events is important because large 
flares, especially on young pre-main sequence stars, can pro- 
vide info rmation about the e xte nt of X-ray coronae (s ee for 
example. [Mullan et al1l2006l and lGetman et al.ll2008bh . Ion- 
ization by large X-ray flares can significantly influence the 
chemistry and turbulence (via the magneto-rotational insta- 
bility) of circumstellar discs, which can have profound effects 
on accretion, dust se ttling, protoplanet mig r ation and other 
phys ical processes l)llgner fe Nelsonl I2006L iFeigelson et al.l 
I2010T I 

Several interpretations of multiple peaked flares, often 
bas ed on solar analog ies, have been proposed. For exam- 
ple, iReale et alJ l|2004h observed an X-ray flare on Proxima 
Centauri that showed two distinct peaks in its lightcurve. 
They concluded that the second peak was probably pro- 
duced thxmudijijrimn^ second loop system. Sim- 
ilarly, |L62ez^SantiagoIet^al] (|2010l ) reported the observation, 
by XMM-Newton, of the unusually long (~36ks) rise phase 
of a flare on a young M star in the TW Hya association. 
They interpreted this rise phase as being a result of the su- 
perposition of multiple flares in separate loop systems. This 



interpretation of stellar flares with similar morphologies is 
common in the literature ( see for example IPillitteri et al.l 
120051 . IPandev fc Shighll2008r i . 



In this work, we consider an interesting geometric al- 
ternative to the explanations given above. In this alterna- 
tive, atypical flare morphologies are not the result of mul- 
tiple flare events, but are the result of the time variable 
eclipsing of the flaring coronal plasma caused by the rota- 
tion of the host star. Previous studies have used this inter- 
pretation to explain the morphologies of stellar flares us- 
ing eclipsing by th e flares' host stars ( Ski nner et al.l [l997l . 
IStelzer et al. I I1999T) or by a companion star in eclipsing bi- 
nary systems (Schmitt &: Favatalll999l, Schmitt et al.ll2003l . 
ISanz-Forcada et al.ll2006l . ISanz-Forcada et al.ll2007l ). In this 
work, we explore the eclipsing interpretation within the 
framework of a single loop model. Based on the solar anal- 
ogy, it has recently been argued that it is unlikely that the 
large stellar flares considered here tak e place within a sin- 
gle magnetic loop (jGetman et al.ll201l[ h However, the single 
loop assumption is often taken as a good approximation in 
situations where there is a single dominant loop within a 
complex loop system. 

The way in which eclipsing can produce flares with atyp- 
ical morphologies can be seen in the following hypothetical 
situation, shown in Fig.[T] in which a flare appears to show a 
double peaked morphology. In this example, the beginning of 
the rise phase (flare phase equal to 0.0) occurs when the flar- 
ing magnetic loop is on the limb of the stellar disc. Initially, 
as chromospheric plasma is evaporated into the corona, the 
visible emission measure increases. However, as the star ro- 
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tates, flaring plasma is rotated out of view, resulting in a 
shallower rise from phases 0.0 to 0.25. As the rate at which 
flaring plasma is eclipsed becomes equal to and then exceeds 
the rate at which plasma is added to the corona, an initial 
peak is seen (flare phase equal to 0.23) followed by a grad- 
ual decay in the visible emission measure. However, the host 
star is at an inclination angle such that the flaring loop is 
never totally eclipsed. The flare's rise phase ends at flare 
phase equal to 0.3 and the decay phase begins. This, how- 
ever, is not seen in the visible emission measure curve. As 
the eclipsed section of the flaring magnetic loop begins to ro- 
tate back into view, a second increase in the visible emission 
measure is seen. As the rate at which the flare's total emis- 
sion measure decreases, it equals and then exceeds the rate 
at which eclipsed plasma is rotated back into view. Thus, a 
second peak followed by a second decay phase is seen in the 
flare's emission measure curve. 

Although individual stellar flares have been studied 
in detail, in the last few years it has become possible to 
study large homogeneous samples of flares. The largest 
such study is the Ch andra Orion-Ultradeep Project (COUP; 
iGetman et"alll2005h . In 2003, the Chandra X-Ray Observa- 
tory provided 13 days of near-continuous observations of the 
member s of the Orion Nebu lar Cluster. Using these obser- 
vations, [Getman~etJil] (|2005l ) identified 1616 X-ray sources, 
of which ~1400 sources were confirmed as members of the 
Orion star forming region and the majority of the rest be- 
ing backgroun d quasars seen thro ugh the molecular cloud 
of the region JCetman et al.l l2005l. Using the COUP data, 
IGetman et alj j2008ah reported the detection of 216 bright 
flares on 161 of these stars using the condition that a 'bright 
flare' is any event that has a peak count rate above four times 
the characteristic (quiescent) count rate for the host star. 
This sample consists of the longest, brighte st and hottest 
flares d etected during the COUP mission. IGetman et al.l 
(2008aj) derived emission measures, flare dur ations, flare 
loop l engths (using the hydrodynamic models of lReale et al.l 
|l993)) and other parameters which will be used exten- 
sively in this paper. A scheme was defined that classified 
flares based on their lightcurve morphologies. Each of the 
216 flares were classified as typical (84), double (8), step 
(38) , slow-rise top-flat (2 0), other (24) or incomplete (42) 
(see IGetman et alJ l2008al for precise definitions and exam- 
ple flare morphologies). 

In this paper, we consider the eclipsing interpretation 
of atypical stellar flare morphologies. More specifically, we 
ask whether such morphologies can be produced through the 
eclipsing of typical flares and then ask to what extent the 
atypical COUP flares are likely to have been produced in 
this way. In Section 2, we describe the simple flare model 
which is used throughout this paper. In Fig. we use three 
examples of COUP flares to show that a range of atypical 
flare morphologies can be explained by eclipsing. In Section 
we consider the distribution of atypical COUP morpholo- 
gies and compare it to a similar modelled set of flares. In 
Section [4] we consider the effect that eclipsing of flares can 
have on the determination of flare loop lengths. Finally, in 
Section [5] the main results and our conclusions from the 
paper are summarised. 




Figure 2. Example of an elliptical loop used to illustrate the 
parameters that define a flare's geometry. H is the loop height, 
W is the loop width, 7 is the angle between the plane in which 
the loop is contained and the stellar rotation axis, and 9q and 4>0 
are the latitude and longitude of the centre of the loop on the 
stellar surface. 



2 FLARE MODEL 

Our flare model involves the following assumptions 

• The flare's emission is the result of a single event con- 
sisting of a rapid rise phase followed by a slow exponential 
decay of the flare's emission measure with time. When mod- 
elling the variation of emission measures with time, we do 
not consider the physical mechanisms that are responsible 
for triggering the flare. 

• The flaring plasma is completely contained within a sin- 
gle static magnetic loop with a uniform plasma density along 
its entire length. 

• The geometry of the flaring loop is described by an 
ellipse with its centre located on, and its major axis perpen- 
dicular to, the stellar surface. The thickness of the flaring 
loop is assumed to be a negligible fraction of its length. 

The first assumption, which is made throughout this pa- 
per, means that any deviation from a simple flare lightcurve 
morphology in the modelled flares can only be a result of 
eclipsing of the flaring plasma. 

The flaring loop geometry is thus characterized by the 
following five quantities: the height of the apex of the loop 
(i.e. the semi-major axis of the ellipse) (H), the width of 
the loop (i.e. the semi-minor axis of the ellipse) (W), the 
latitude and longitude of the centre of the ellipse (80, (j>o), 
and the angle between the plane of the ellipse and the star's 
rotation axis (7). These quantities can be seen in Fig. [2] 
It is worth emphasising that the term 'width' in this case 
refers to the length of the semi-minor axis of the ellipse and 
not the more common definition of the length between the 
two loop footpoints along the segment of the great circle that 
connects them. For the purposes of this paper, the difference 
between these two definitions is not important. The other 
parameters that can determine the effects of eclipsing are 
the stellar inclination angle, the stellar rotation period, the 
flare's duration and peak emission measure. 
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Figure 3. Emission measure vs time curves for the observed flares on COUP sources 66 (left), 649 (middle) and 942 (right) which represen t 
good examples of atypical COUP flares. The asterisks show the observed COUP emission measure data, given bv lGetman et al.l {2008a). 
The shaded area represents the times when some portion of the magnetic structure containing the flare was eclipsed by the host star. 
The solid lines show the intrinsic emission measure curves, EMtotif), which have been fitted to the observed emission measure curves 
in the region outside the shaded area; this represents what the flare may have looked like had the flare always remained in view. The 
dashed line shows an eclipsed version of the same modelled flare which gives the best fit to the observed flare's emission measure curve. 



Under the assumption that the plasma density is uni- 
form over the length of the flaring loop, the visible emission 
measure, EM V i S (t), can be expressed as 

EM ms {t) = EM tot (t) (1) 

\ Vtot J 

where EM tot (t) is the flare emission measure curve that will 
be seen if the entire flaring loop is visible throughout the 
duration of the flare and V V i 3 (t) /Vtot is the fraction of the 
volume of the flaring plasma that is visible at any given 
time t. The quantity V v is(t) /Vtot is calculated at each time 
t by considering a series of points equally spaced along the 
length of the flaring loop. Under the assumption that the 
loop thickness is small, the fraction of the flaring loop vol- 
ume that is visible at this time is approximately equal to 
the fraction of these points that are visible. We give details 
of how to determine whether a point on a flaring loop is 
eclipsed or visible in Appendix [X] The methods used for 
chosing the function EMtot(t) is described in Paragraph 1 
of Section \3. II 

The most obvious source of eclipsing of stellar flares 
comes from the host stars which we assume to be opaque 
spheres. However, other sources of eclipsing may be present. 
In pre-main sequence stars, these may be binary compan- 
ions, circumstellar discs, accretion columns extending from 
a circumstellar disc to the stellar surface, planets at small 
radii and flare associated prominences. In this paper, the 
only sources of eclipsing that we consider are host stars, 
circumstellar discs and flare associated prominences. 

We model circumstellar discs as opaque discs with 
smooth inner edges located at the equatorial corotation radii 
(Rao = (GM./W 2 ) 1 , where cj is the angular velocity of 
the stellar surface at the equator) of their host stars. The 
modelled discs are assumed to be flat and to lie in the equa- 
torial plane. The possibility of more complex discs is not 
considered here although it should be noted that a warped 
circumstellar disc could have a significant effect on a flare's 
lightcurve, particularly if the stellar inclination is such that 
a warped inner disc periodic ally obscures the view to the 
star (e.g. lAlencar et al.|[201ol ). 



We model prominences as opaque spheres that sit above 
the apex of flaring loops. Thus, a prominence is characterised 
by its height above the flaring loop and its radius. We take 
all prominences to be spheres of radius 0.57?*, the centres 
of which have heights above the apex of the flaring loops of 
0.55i?». The prominences thus cover 25% of the stellar disc, 
which is similar to the estimated pr ojected areas of promi- 
nences on AB Dor and Spe edy Mic (jCollier Cameron et al.l 
1 1990) ; iDunstone etal]|2006h . 



3 THE NUMBER OF ECLIPSED FLARES 

In this section, we analyse the entire COUP sample in order 
to determine how many flares have been eclipsed. We define 
non-eclipse candidate flares as flares that show a single rise 
followed by a single decay in their emission measure vs time 
curves (this includes bot h the 'typical' and th e 'slow-rise 
top-flat' flares defined bv lGetman et ai1l2008al ). We define 
eclipse candidate flares as those that display sudden short 
duration decreases followed by increases in their emission 
measures. 

An eclipse that has a duration comparable to 
the duration of the flare will generally only result 
in a less luminous flare without a noticeably atypi- 
cal morphology or a flare that is not visible at all. 
For this reason, we expect that such dips should 
be found predominantly on long duration flares and 
rapidly rotating stars. 

In order to illustrate the effects that eclipsing 
can have on typical flares, we show in Fig. [3J three 
examples of eclipse candidate flares from the COUP 
sample. We demonstrate that an eclipse candidate 
flare can be produced by the eclipsing of a typical 
flare (i.e. flares consisting of a single rise phase fol- 
lowed by a slow exponential decay) by fitting model 
eclipsed flares to these observed emission measure vs 
time curves (for details, see the caption of Fig. [3|). 
We note, however, that while the forward problem of 
varying the model flare parameters to fit the obser- 
vations is quite straightforward, the inverse problem 
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Figure 4. Flare on COUP source 871 which our standard 'typical' 
emission measure curve is based on. The full line represents the 
best fit to this curve which is used as the standard flare emission 
measure vs time curve. 



of recovering the true flare parameters solely from 
the observations is in general not possi ble. 

Of the 216 COUP flares catalogued by iGetman et al.l 
(2008a|), we identify 62 (29%) eclipse candidates. This num- 
ber, given different levels of scepticism by the examiner, may 
be between 31 (14%) and 71 (33%). For the rest of this pa- 
per, we will take the value of 62 as the number of eclipse 
candidate flares in the COUP sample. 

Table [T] gives average values for several COUP flare 
and host star paramet ers, derived from parameters given by 
IGetman et al.l (|2008al ) , for eclipse candidate and non-eclipse 
candidate flares separately. We have attempted to estimate 
by eye any decreases in the flare durations derived from the 
visible emission measure curves that might have been caused 
by eclipsing. This is only possible for flare emission measure 
curves which have been broken into two parts by large tem- 
porary eclipses. It can be seen that contrary to expectations, 
the average flare durations as a fraction of host star rotation 
period is shorter for eclip se candidate fla r es whe n the orig- 
inal flare durations from IGetman et al.l (|2008a ) are used. 
When the durations are calculated assuming eclipsing has 
occurred, this is no longer a problem because the durations 
are always significantly longer than their original values. It 
can also be seen that the average peak emission measures are 
lower for eclipse candidate flares than for non-eclipse candi- 
date flares which is consistent with eclipsing hypothesis. 



3.1 Modelling the set of COUP flares 

Given that we classified 62 of the COUP flares as eclipse 
candidates, we now use our simple flare model to calculate 
the number of eclipse candidates that we would expect to 
see in the COUP flare sample. We assume that all flares are 
produced by the same energy release process acting within 
our simplified loop geometry, and therefore have similar typ- 
ical intrinsic emission measure vs time curves (i.e. EM to t(t) 
in Eqn.Q}. We therefore choose the flare observed on COUP 
source 871 from the COUP sample (shown in Fig. |4| as the 
intrinsic flare profile. We then use this profile to produce 
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Figure 5. Histograms showing the distributions of relevant stellar 
and flare parameter s deriv ed from the data given in Tables 1- 
3 of IGetman et al.l d2008all . The histograms show, from top to 
bottom, stellar rotation periods, stellar masses, stellar radii, host 
star characteristic emission measures, flare durations, flare loop 
lengths and flare peak emission measures. Dashed lines show the 
locations of the mean values. 
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many simulated flares by scaling the peak emission measure 
and the flare duration in the following way 

EM tot , alm (t) = EM tot ,sn {t acaled ) x f EMp^l) ( 2 ) 

\ fiMmai,871 J 

where 

''scaled = t X — yo) 
E871 

where the time t is zero at the beginning of the impulsive 
phase of the flare, EMtot,87i, EM m ax,s7i and ts7i are re- 
spectively the intrinsic emission measure vs time curve, the 
maximum emission measure, and the duration for our stan- 
dard flare and EM to t,sim, EM m ax,sim and t sim are similar 
quantities for the simulated flare. We choose the values of 
EMmax.sim and tsim randomly using the model described 
below in Section 13.21 We note that although here we have 
focused on the emission measure vs time curve of the flare 
observed on COUP source 871 as our standard typical flare, 
repetitions of our analysis as discussed below using a differ- 
ent typical COUP flares yielded no significant difference in 
our results. 

In order to determine the visible emission measures 
curves (EM v i a (t)) for the simulated flares we also need to 
specify the nine parameters that determine the loop geom- 
etry and position and the stellar rotation rate. Each of the 
randomly chosen parameters discussed above is chosen based 
on probability distributions that best approximate the dis- 
tributions of these parameters in the COUP sample. Where 
it is not possible to use observed distributions of parameters, 
reasonable assumptions, discussed below, must be made. 
With these parameters we calculate the fraction of the flar- 
ing loops that are visible as a function of time (i.e. the frac- 
tion V V is/Vtot from Eqn. [T]) and using Eqn. [TJ we calculate 
the visible emission measure curves for each flare. 

It is importa nt in these calculation s to define flares in 
the same way as iGetman et "all (|2008aT l in order that the 
results can be reliably compared. For this reason, we define 
a flare as any energetic event with a peak emission measure 
exceeding four times the characteristic emission measure of 
the host star. In this way, flares that have been eclipsed to a 
an extent that they would not have been classified as bright 
flares in the COUP sample are discarded. 

3.2 Probability distributions for flare parameters 

In Section 2, we listed the nine geometric and temporal pa- 
rameters that can affect, through eclipsing, the soft X-ray 
lightcurve morphologies of stellar flares. Another factor con- 
sidered in this paper is the existence of other opaque mate- 
rial that can act as alternative sources of eclipsing. The two 
other sources of eclipsing considered here are circumstellar 
discs and flare associated prominences. In order to calcu- 
late the radii of the inner edges of the circumstellar discs, 
which we assume to be at the equatorial corotation radius, 
we must also model the stellar masses and radii. Thus, for 
the purposes of this paper, we must model eleven probability 
distributions. 

Fig. [5] shows histograms for seven of the p a ramete rs de- 
rived using the data given by I Get man et al.l j2008a). For 
these distributions, we ignore data from eclipse candidate 
flares. With the exception of the stellar rotation periods, all 



these parameters can be modelled using log-normal distri- 
butions parametrised by their mean, /x, and variance, a 2 . 
The means and variances for these six parameters are given 
in Table [2] We assume that the stellar rotation periods have 
values that are evenly distributed between 0.1 and 11 days. 

The other parameters that need to be estimated in or- 
der to model the ONC flares are the starting longitudes of 
the flaring loops, <j>o, the orientation of the loops, 7, the in- 
clination angles of the stellar rotation axes to an observer's 
line-of-sight, 6 V i £W , and the colatitudes of the centres of the 
flaring loops, 8i oop . The former two are taken to have val- 
ues that are evenly distributed over all possible values. The 
latter two are taken to have a higher probability for values 
near the equator based on the probability density function 
pdf(8) — i cos8i at , where 8i at is the latitude. 

In these calculations, the presence of circ umstellar discs 
around some of the stars is also considered. IGetman et alj 
(2008a) derived near infrared colour excess (A(_ff — K s )) 
values for 140 of the flare host stars and used the condition 
A(H — K s ) > —0.06 mag as a good indicator for the pres- 
ence of circumstellar discs. Of these 140 stars, 53 indicate 
the presence of a circumstellar disc. Thus, in the flares sets 
considered in the next section, where circumstellar discs are 
considered, each flare has a probability of 0.38 of having 
occurred on a star that has a disc. 

It is also necessary to calculate the heights (H) and 
widths (W) of flaring loops when the only available informa- 
tion are the loop lengths. For this reason, it is then necessary 
to assume a plausible relation between the heights and the 
widths of flaring loops. We assume the relation that would 
be expected for a potential arcade with a maximum width 
of Wm ax- This relation can be found in lBrowning fc Priest! 
( 1986) and is given by 

° xp (-£) =cos (2i^) (4) 

where W max is taken to be equal to 0.9R* (if the value for 
Wmax is larger than R* the largest flare loops would not 
touch the stellar surface). It is important to point out that 
even though we use the height-width relation for magnetic 
loops in a potential arcade, throughout this paper, the actual 
loop geometries are ellipses. 

3.3 Results 

Using the method described above, we model a set of 10,000 
flares. In this section, we analyse these flares using the same 
method of visual inspection to select eclipse candidate flares 
as was used to analyse the COUP sample. In the initial re- 
sults presented here, circumstellar discs and flare-associated 
prominences are not considered. These are included in sep- 
arate results presented at the end of this section. 

Due to eclipsing, not every modelled flare has a peak 
emission measure that is large enough to be classified as a 
bright flare in the COUP observations. In order to obtain 
10,000 flares which make it through the observa tional selec- 
tion criteria adopted by IGetman et al.l (l2008aT ) it is neces- 
sary to model 10,878 flares in total. By inspecting the entire 
sample of modelled flares, we classify only 7.0% as eclipse 
candidates. This small value is to be contrasted with the 
larger number of 29% eclipse candidate flares in the COUP 
sample. 
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N tot 
(1) 


Prot 

(days) 
(2) 


EM pk 
(10 53 cm" 3 ) 
(3) 


L 

(10 10 em) 
(4) 


N disc /Ntot 
(5) 


*/ 

(days) 
(6) 


(COUP) 
(tf/Prot) 

(7) 


tf (Modified) 
(days) (t f /Prot) 
(8) (9) 


Eclipse Candidates 


62 


7.08 


33.92 


42.71 


0.33 


1.28 


0.27 


2.36 


0.47 


Non-Eclipse Candidates 


154 


6.57 


46.87 


29.22 


0.44 


1.59 


0.36 


1.59 


0.36 


All 


216 


6.71 


43.15 


33.09 


0.40 


1.50 


0.34 


1.85 


0.41 



Table 1. Average values of flare and stellar parameters in the COUP set reported by iGetman et all d2008al) . The data are presented 
for the entire set (bottom row) and separately for the flares that have been classified as eclipse candidates (top row) and non-eclipse 
candidates (middle row). The columns correspond to: Col. (1): the number of flares in each category. Col. (2): the host star rotation 
periods. Col. (3): the peak emission measures. Col. (4): the derived loop lengths. Col. (5): the fraction of flares on st ars which show 
evidence that they possess a circumstellar disc. Cols. (6,7): the visible flare durations reported bv lGetman et aL (2008a) as an absolute 
value and as a fraction of their host star's rotation period. Cols. (8,9): the flare durations after we attempt to correct for decreases in 
their visible durations due to noticeable eclipses. 



log w (tf)(days) 


log 10 (L)(R t ) 


log 10 (M*/M Q ) 


log 10 (EM pk )(cm 3 ) 


logio(EM char )(cm 3 ) 


log 10 {R t )(10 10 cm) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


fi 0.037 


0.30 


-0.07 


54.37 


53.71 


1.25 


cr 2 0.13 


0.24 


0.13 


0.26 


0.25 


0.05 



Table 2. Mean and variance values for the six flare and host star parameters that are to be approximated as log-normal distributions. 
The columns correspond to: Col. (1): flare durations. Col. (2): Flare loop lengths. Col. (3): Mass of host stars. Col. (4): Flare peak 
emission measures. Col . (5): H ost star characteristic emission measures. Col. (6): Host star radii. The values are calculated using data 
given by Getm an et al . (2008a). which are plotted as histograms in Fig. [5] 




Flare Phase Flare Phase Flare Phase 

(d) (e) (f) 

Figure 6. Six examples of modelled flares produced through eclipsing of the flaring plasma where the solid lines show the flare emission 
measure curves prior to eclipsing (i.e. the emission measure curve that would have been observed had the event always remained in view) 
and the dashed lines show the visible (i.e. eclipsed) emission measure curves. The latter arises by allowing the flare to enter or exit from 
rotational eclipse. These examples show that it is likely that atypical flare morphologies will be produced at random given a large sample 
of flares. 
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Ntot 


Prat 


EM p k (visible) 


EMpk (original) 




(visible) 




(original) 






(days) 


(10 53 cm,- 3 ) 


(10 53 cm,- 3 ) 


(days) 


(tf/Prot) 


(days) 


(tf/Prot) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Eclipse Candidates 


698 


2.07 


49.0 


59.7 


1.51 


1.91 


2.26 


2.63 


Non-Eclipse Candidates 


9302 


5.77 


46.5 


49.6 


1.31 


0.47 


1.33 


0.48 


All 


10,000 


5.55 


46.6 


50.1 


1.32 


0.57 


1.39 


0.61 



Table 3. Average values of flare and stellar parameters in the set of modelled flares. The data are presented for the entire set (bottom 
row) and separately for eclipse candidates (top row) and non-eclipse candidates (middle row). The columns correspond to: Col. (1): the 
number of flares in each category. Col. (2): the stellar rotation periods. Col. (3): the visible peak emission measures. Col. (4): the actual 
peak emission measure. Cols. (5,6): the visible flare durations as an absolute value and as a fraction of the host stars' rotation periods. 
Cols. (7,8): the actual flare durations prior to eclipsing. 



Six examples of flare emission measure curves which 
have been affected by eclipsing are shown in Fig. [6] We de- 
termine that eclipse candidates are likely to occur in a sam- 
ple of 216 flares under the conditions present in the ONC. 
Thus, we conclude that it is likely that a number (some, 
although not necessarily all) of the eclipse candidate flares 
in the COUP sample have been produced by the rotational 
eclipsing of typical flares. However, it is unlikely that the 
entire sample has been produced in this way. Given that 
the probability of one flare being an eclipse candidate is 
0.07, using the Binomial distribution, we estimate that the 
probability of 62 flares being eclipse candidates in a sam- 
ple of 216 flares is approximately 10 -25 . In order to explain 
the large number of atypical COUP flares, it is thus nec- 
essary to assume other physical mechanisms, such as mul- 
tiple heating events in a single flaring loop, multiple flares 
whose lightcurves have been superimposed, or stellar ana- 
logues of solar coronal arcades where a reconnection event 
triggers subsequent events and associated flares/heating of 
neighbouring loops. 

The eclipse candidates in the modelled set of flares, how- 
ever, do not represent the full sample that have undergone 
eclipsing. The fifth and sixth flares shown in Fig. [5] show ex- 
amples of flares that have been eclipsed but still have 'typ- 
ical' emission measure curves. A total of 63% of the mod- 
elled flares have been partially eclipsed for at least a frac- 
tion of their durations and 49% of the flares have their peak 
emission measures reduced. However, in most cases, eclips- 
ing leads to an insignificant decrease in the visible emission 
measure value. This can be seen in Fig. [7] Of the modelled 
flares, 6% show significant decreases in their visible dura- 
tions. Therefore, we expect that analysis of such flares may 
lead to derived flare parameters that are different from the 
true physical properties of the magnetic structure containing 
the flare; we explore this point further in the next section. 

One type of flare morphology seen in the COUP sample 
but not in the fla res modelled h e re are the slow-rise top-flat 
flares, defined bv lGetman et all {2008a). Such morphologies 
can be produced through eclipsing, as can be seen in the 
second peak in the fourth example shown in Fig. [6] However, 
this peak has not made it through the selection criterion for 
flares so it is not counted in the sample of modelled flares. 

In TableO we give average values for flare and host star 
parameters, both before and after eclipsing has been taken 
into account, for eclipse candidate and non-eclipse candidate 
flares separately. It can be seen that eclipsing causes a reduc- 



tion in the average durations and peak emission measures. It 
can also be seen that longer duration flares on faster rotating 
stars are more likely to be eclipse candidates, as expected 
for flares randomly distributed in latitude and longitude. 

One reason why there may be more eclipse candidates 
seen in the observed COUP sample than in the modelled 
flare set could be that the COUP flares are being eclipsed 
by circumstellar discs or flare-associated prominences. To 
investigate to what extent this may be the case, we repeat 
the above calculations to produce two more sets of 10,000 
flares. In the first set, we assume that circumstellar discs 
are present ar ound 38% of the host s tars; the same fraction 
as determined iGetman et all (|2008bl ) based on Spitzer H-K 
excess emission. With this assumption, we find that 6.4% are 
eclipse candidates. In order to produce 10,000 visible flares, 
it was necessary to produce 11,689 flares in total. In the 
second set, we assume that prominences are present above 
the apex of each flaring loop. With this assumption, we find 
that 7.7% are eclipse candidates. In order to produce 10,000 
visible flares, it was necessary to produce 11,612 flares in 
total. Thus, it is clear that even with circumstellar discs 
and flare-associated prominences we are not able to explain 
all of the eclipse candidates seen in the COUP sample with 
eclipsing of single magnetic loops only. 



4 FLARE LOOP LENGTHS 

A common method for determining the loop length of an 
unresolved stellar flare involves the comparison of the flare's 
emission measure and temperature data with hydrodynami c 
flare models. This method, detailed by iReale et al.1 (|l997l ). 
defines the loop half-length as 



L(cm) 



T d (ks)y/T pk (MK) 
3.7 x W- 4 F(Q 



where 



F(0 



0.63 
C - 0.32 



+ 1.41 



r pk = 0.068T o ^ 



(5) 

(6) 
(7) 



where Td(ks) is the time that it takes for the flare's emission 
measure to decay by a factor of e (the e- folding timescale), 
Tpfe(MK) is the temperature at the apex of the flaring loop 
when the flare's emission measure is at its peak, T (, S (MK) 
is the observed average loop temperature at this time and £ 
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Figure 7. Histogram showing the visible peak emission measures, 
EMpk, as a fraction of their uneclipsed values, EM v ^ tot , for the 
set of modelled flares. Only the 105 flares which had their peak 
emission measures decreased by eclipsing are included here. 



is the gr adient of the decay phase of the log(EM?) - logT 
plot. See lReale et al.l (|l997l ) for full details (for a discussion 
on the validity of th e single loop model, see Appendix A of 
iGetman et ail (l201rT )l. 

In order to explore the effect that eclipsing can have on 
the loop lengths derived from this method, which has been 
com monly employed in the analysis of flares on young stars 
fe.g. lFavata et ail (|2005t )l we take the emission measure and 
temperature data for the 'typical' flare seen on COUP source 
871 and produce a large set of 10,000 flares using the method 
described in Section T3.2I In these calculations, however, we 
only pick at random the stellar parameters (e.g. radius, ro- 
tation period) and the locations and orientations of the flar- 
ing loops. The parameters specific to the flare (e.g. peak 
emission measures, lo op length, flare duratio n) are kept at 
the values derived by IGetman et all |2008al ). We then ap- 
ply the loop length analysis to the uneclipsed and eclipsed 
flare emission measure curves, assuming that eclipsing has 
no effect on the determined temperatures. 

The effect of eclipsing is shown in Fig. [8] The loop 
length based on the uneclipsed emission measure curve is 
14.9 x 10 10 cm. After eclipsing has been taken into account 
for a large set of flares, the average calculated loop length 
for that set is increased to 16.8 x 10 10 cm. This simple esti- 
mate suggests that in most cases eclipsing has little effect 
on derived loop lengths. 



5 SUMMARY AND CONCLUSIONS 

Although most stellar flares have typical soft X-ray 
lightcurve morphologies (i.e. a single rapid rise followed 
by a slow exponential decay), many flares have atypical 
morphologies. Many of th ese show multiple pea ks or small 
dips in their lightcurves (|Getman et al,ll2008al l Based on 
solar analogies, such flares are often interpreted as being 
the result of multiple heating events in the same flaring 
loop or the super position of separ a te overlapping flares 
(for example, see iReale et al.l |2004|; iPillitteri et al.l 120051 ; 



iPandev fc Singbl 120081 ; iLopez-Santiago et al.ll2010h . In this 
paper, we have considered an alternative geometric inter- 
pretation in which these atypical flare morphologies are pro- 
duced by the eclipsing of flaring plasma due to the rotation 
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Figure 8. Histogram showing the derived lengths of flaring loops 
after eclipsing has been taken into account (flares which have 
not undergone any eclipsing are included). The solid vertical line 
shows the length that would have been derived for all of the flares 
if no eclipsing had occurred. The dashed vertical line shows the 
log of the average derived loop lengths (not the average of the 
log). 



of the host star. This interpretation has been considered 
for individua l flares in previous studies ll Skinner et al.lll997l 



i Stelzer et all I1999I, ISchmitt fc Favatal Il999l. ISchmitt et all 
booi ISanz-Forcada et alll200rj| . ISanz-Forcada et al.ll2007f l . 

Using dat a from t he Chandra Orion Ultra-deep Project, 
IGetman et all (|2008al ) identified 216 stellar flares on 161 
PMS stars. As the COUP sample contains a range of flare 
morphologies, in this paper, we have used it to explore the 
eclipsing interpretation. We analysed the entire COUP sam- 
ple by eye to determine which of them are eclipse candidate 
flares. In Fig. [3] we took three examples of these and showed 
that their emission measure vs time curves can easily be 
produced by the eclipsing of typical flares. In Section [3] we 
showed that the entire COUP sample contained 62 (29%) 
eclipse candidates. However, by producing a large modelled 
set of flares similar to the COUP sample, we showed that 
although 63% of the modelled flares underwent eclipsing, 
this was detectable in only 7.0% of them. In Section [3] we 
showed that eclipsing can effect the derived loop lengths for 
flares, but in most cases such an effect is negligible. 
Our conclusions from these results are as follows 

• The time variable eclipsing of stellar flares contained 
within single magnetic loop structures can produce the atyp- 
ical morphologies observed in the COUP sample. Thus, 
given a flare with an atypical lightcurve, it is not necessary 
to invoke unusual physical mechanisms to account for the 
flare's morphology. However, it should be noted that eclips- 
ing is much more likely to cause an atypical morphology on 
longer duration flares and on more rapidly rotating stars. 

• However, the observed frequency of eclipse candidate 
flares in the COUP sample is far higher than we would ex- 
pect if eclipsing was the only mechanism by which atypical 
flares were being produced. Thus, alternative physical mech- 
anisms, such as the stellar analogies of solar coronal arcades, 
must be responsible for most of the atypical COUP flares. 

• Even in cases where an observed flare has a typical mor- 
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Observer 
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Figure Al. The various coordinate systems used to determine 
whether or not an object has been eclipsed. In this image, the y 
and y' axes are going into the page. 

phology, it is not possible to know from the flare lightcurve 
alone whether or not eclipsing has taken place. 

• Eclipsing is unlikely to have a significant effect on de- 
rived loop lengths. 
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APPENDIX A: CONDITIONS FOR ECLIPSING 

In this appendix, we describe the method used in this paper 
to determine whether or not a point is being eclipsed by 
either the host star, a circumstellar disc or a prominence. 

It is easiest to set up and evolve a stellar system in a 
spherical polar coordinate system (r,9i a t,<j))- However, the 
conditions for eclipsing are simplest when the system is rep- 
resented in Cartesian coordinates with either one axis point- 
ing along the line of zero longitude and latitude (x,y,z) or 
pointing towards the observer (x' , y' , z'). These three coordi- 
nate systems can be seen in Fig. I All In all three coordinate 
systems, the origin is at the centre of the host star. The 
transformations between these three coordinate systems are 
as follows. 



x — r COS (j) COS 8i at 


(Al) 


y — r sin cj> cos 9 ta t 


(A2) 


z = rsmOiat 


(A3) 


— z sin 9iat,o + x cos 9i a t,o 


(A4) 


y =y 


(A5) 


= z cos Biat,o - x sin 0i a t,o 


(A6) 



where 9i at is the latitude and defined such that it has 
values between —90° and 90° and the observer is located 
at (r,9iat,4>) equal to (oo, Oi a t,o, 0) and (x',y',z') equal to 
(OO.0.0). 

Consider a point which does not lie within the star, 
the disc, or the prominence, and has coordinates (r, 9i a t, 0), 



(x,y,z) and (x',y',z') in the different coordinate systems. 
This point is only visible if it is not eclipsed in all three of 
the following conditions. 

Condition 1: The point is eclipsed by the host star with 
radius ii* if 

x < (A7) 
y' 2 + z' 2 < R 2 (A8) 

The latter condition is only met if the point coincides 
with the disc of the host star in the plane of the sky and 
the former condition tests if the point is behind or in front 
of the host star. 

Condition 2: The point is eclipsed by a circumstellar 
disc with an inner hole with a radius of Rtrunc, which we 
assume to be the equatorial corotation radius, if 

9iat x 9 lat , < (A9) 

sin 2 9iat,oz' 2 + y' 2 > R 2 runc (A10) 

As 9iat and 9i a t,o are defined such that they are positive 
in one hemisphere and negative the the other, the first of 
these condition is only met if the point and the observer are 
in opposite hemispheres of the star. The second condition 
then tests if the point is visible through the inner hole of 
the disc as projected onto the plane of the sky. 

Condition 3: The point is eclipsed by a prominence of 
radius R p with its centre at (r p ,9iat,p,<f>p), (x p ,y p , z p ) and 
(*£p7 y P i z p ) if 

x < x' p (All) 

(y - y'p? + - 4) 2 < Rl (Ai2) 

These conditions are the same as the conditions for 
eclipsing by the host star above with corrections for a sphere 
with a different radius and which is not centred at the origin 
of the coordinate system. 
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